Measurements of J/V' and il>(2S) decays into AAtt" and AArj 
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Using 58 million J/ip and 14 million Tp{2S) events collected by the BESII detector at the BEPC, 
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branching fractions or upper limits for the decays J/ip and tp{2S) — > AAtt^ and AA?7 are measured. 
For the isospin violating decays, the upper limits are determined to be B{J/'^ — » AAtt") < 6.4 x 10~® 
and B{ip{2_S) -» AAtt") < 4.9 x 10"^ at the 90% confidence level. The isospm conserving process 
J/tp — > AAr; is observed for the first time, and its branching fraction is measured to be B{J/tp — > 
AAt?) = (2.62 ± 0.60 ± 0.44) x lO^", where the first error is statistical and the second one is 
systematic. No AAr/ signal is observed in Tp{2S) decays, and B{ip{2S) — > AA?;) < 1.2 x 10~* is set 
at the 90% confidence level. Branching fractions of J/xp decays into E+tt^A and S~7r"'"A are also 
reported, and the sum of these branching fractions is determined to be B{J/%I} E+tt" A + c.c.) = 
(1.52 ±0.08 ±0.16) X 10"^. 

PACS numbers: 13.25. Gv, 12.38.Qk, 14.20.Gk, 14.40.Cs 



I. INTRODUCTION 

Several charmonium decay modes containing AA pairs 
have been reported [1-8]. Among these decays, the 
isospin violating process J/tp ^ AAtt" has been stud- 
ied by DM2 [3] and BESI [4], and its average branch- 
ing fraction is determined to be B{J/tp AAtt^) = 
(2.2 ± 0.6) X lO"'^ [9]. However, the isospin conserving 
process J/tp AAr] has not been reported, and there are 
no measurements for AAtt^ and AAr] decays of ip{2S). 

In this paper, we study J/tp and tp{2S) — > AAtt", AArj 
using 58 M J/tp events and 14 M '0(25') events taken with 
the BESII detector at the BEPC storage ring. We find 
that the J/tp — > AAtt" branching fraction is much smaller 
than those measured by DM2 and BESI. In addition, we 
observe the isospin conserving process J/tp AArj and 
measure its branching fraction for the first time. Analyses 
of AAtt^ and AAt] in tp{2S) decays are also performed, but 
no obvious signals are observed for these two channels. 



II. THE BESII DETECTOR AND MONTE 
CARLO SIMULATION 

BESII is a conventional solenoidal magnet detector 
that is described in detail in Ref. [10]. A 12-layer ver- 
tex chamber (VTC) surrounding the beam pipe provides 
trigger and track information. A forty-layer main drift 
chamber (MDC), located radially outside the VTC, pro- 
vides trajectory and energy loss (dE/dx) information for 
charged tracks over 85% of the total solid angle. The 
momentum resolution is ap/p = 0.0178\/l -I- p"^ {p in 
GeV/c), and the dE/dx resolution for hadron tracks is 
~ 8%. An array of 48 scintillation counters surrounding 
the MDC measures the time-of-flight (TOF) of charged 
tracks with a resolution of ~ 200 ps for hadrons. Radi- 
ally outside the TOF system is a 12 radiation length, lead 
gas-tube barrel shower counter (BSC). This measures the 
energies of electrons and photons over ~ 80% of the total 
solid angle with an energy resolution of as /E = 21%/y/E 
{E in GeV). Outside of the solenoidal coil, which provides 
a 0.4 Tesla magnetic field over the tracking volume, is 
an iron flux return that is instrumented with three dou- 
ble layers of counters that identify muons of momentum 
greater than 0.5 GeV/c. 



In this analysis, a GEANT3 based Monte Carlo (MC) 
simulation program [11] with detailed consideration of 
real detector responses (such as dead electronic channels) 
is used. The consistency between data and Monte Carlo 
has been carefully checked in many high-purity physics 
channels, and the agreement is quite reasonable [12]. 



III. EVENT SELECTION 

The decay channels investigated in this paper are 

J/tP AAf, J/tP AAr/, tP{2S) AAtt^, and 
tp(2S) AAr], where A decays to n~p and 7r° and ry 
to 77. The final states in which we are interested contain 
two photons and four charged tracks (7r+7r~pp). Candi- 
date events are required to satisfy the following common 
selection criteria: 

1. Events must have four good charged tracks with net 
charge zero. A good charged track is a track that 

is well fitted to a helix in the MDC and has a polar 

angle, 6, in the range | cosd\ < 0.8. 

2. The TOF and dE/dx measurements of the charged 
tracks are used to calculate x%id values for the hy- 
potheses that the particle is a pion, kaon, or proton. 
Only the two proton tracks must be identified with 
the requirement that Xpid ^'-'^ ^^"^ proton hypoth- 
esis is less than those for the tt or hypotheses. 

3. Isolated photons are those that have an energy de- 
posit in the BSC greater than 50 MeV and have 
the angle between the photon entering the BSC, 
and the shower development direction in the BSC 
less than 37°. In order to remove the fake pho- 
tons produced by p annihilation and those produced 
by hadronic interactions of tracks with the shower 
counter, the angle between the photon and antipro- 
ton is required to be larger than 25° and those be- 
tween the photon and other charged tracks larger 
than 8°. 

4. The selected events are subjected to four constraint 
(4C) kinematic fits. When there are more than two 
candidate photons in an event, all combinations are 
tried and the combination with the smallest xic 
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retained. The selection requirement on xtc °P' 
timized by maximizing S/y/S + B, where S and B 
arc the expected numbers of signal and background 
events, respectively. 

5. To select A and A, the difference between the 
measured wp mass and the expected mass (M(A)) 
should be less than 10 MeV/c^ (three times the A 
mass resolution). 



IV. EVENT ANALYSIS 
A. J/ip AAtt" 
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Event Selection 



Only events with two good photons are selected, and 
4C kinematic fits under the 777r+7r~pp hypothesis are 
performed. To select clean AA events, we reqiiirc the A 
and A secondary vertices to be reconstructed successfully, 
and the decay lengths of A and A in the x — y plane must 
be larger than 0.05 m. 

Figure 1 shows the AA invariant mass (M(AA)) distri- 
bution after the above selection. The large peak near 2.9 

GeV/c^ is background from J/ip E'^E'', in agreement 
with the expectation from the MC simulation, normal- 
ized to its branching fraction [9], shown as the shaded 
histogram in Fig. 1. To reject such background. M (AA) is 
required to be less than 2.8 GeV/c^. With this selection, 
Fig. 2 shows the xlc distribution for data and Monte 
Carlo simulation. To suppress potential backgrounds, 
xlc < 10 is required. 
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FIG. 1: Distribution of M(AA) for J/ip — » AA77 candidates. 
Dots witli error bars are data, the shaded histogram is back- 
ground from MC simulated J/ip E^E", normalized accord- 
ing to the branching fraction in the PDG, and the dashed 
histogram is the MC simulated J/V' — > AAtt'^ signal, normal- 
ized according to the branching fraction in the PDG. 



FIG. 2: Distribution of xlc for J/tp — » AA77 candidate events 
(solid histogram) and Monte Carlo simulated J/ip AAtt" 
events (dashed histogram). Here, M(AA) is required to be 
less than 2.8 GeV/c^ 
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FIG. 3: Scatter plot of M(7r"'"p) versus M{-w~p) invariant mass 
for (a) J/ip — »■ E"'"7r~A (-|-c.c.) candidate events and (b) MC 
simulation. 



2. Background Analysis 



To explore other possible backgrounds, we generate 
MC events for the following channels: J/ip 7AA, 
S^EO, E(1385)"E(1385)0, E°E°, 2(1530)05°, I:°tt°A + 
c.c, and E+tt^ A + c.c. Only the last two channels give 
significant contributions to the 7r° signal. In particular, 
the decay mode J/ip — > Ti^w^A + c.c, which contains 
AAtt^ with an additional photon in the final state, could 
contribute to the observed number of AAtt^ candidates. 
Because direct measurements of J/ip E°7r°A -|- c.c. 
are difficult, we measure the branching fractions of their 
isospin partners and estimate their branching fractions 
by assuming isospin symmetry. To estimate the contam- 
ination from J/ip ^ E°7r°A + c.c, a high precision mea- 
surement of J/ip S+7r~A -|- c.c. is very important. 
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FIG. 4: Scatter plot of M(77) versus M(7r+p) for (a) J/^) 
E"''7r~A candidate events and (b) MC simulation of J/ip — > 
E^7r~A, both satisfying xlc <15. The central box in the fig- 
ure is the signal region defined by |M(77) — M{tt'^)\ < 0.03 
GeV/c^ and |M(7r+p) - M(A)| < 0.006 GeV/c^ The tt" 
sideband is defined by |M(77) - 0.06| < 0.03 MeV/c^ and 
|M(77) - 0.21| < 0.03 GeV/c^ (The two boxes located above 
and below the signal region), and the A sideband region by 
\M{-K+p) - 1.101)1 < 0.006 GeV/c^ and \M{tv+p) - 1.131)| < 
0.006 GeV/c^ (The two boxes on the left and right of the sig- 
nal region) . The four boxes at the corners are used to estimate 
the phase space contribution. 



30 McV/c2 and \M{t:+p) - M(A)| < 6 MeV/c^), and a 
clear S+ signal is observed. The dashed histogram is 
the background coming from sidebands of 7r° and A. To 
obtain the number of E"*" events, we fit the S+ signal 
with a histogram of the signal shape from MC simulation 
plus the background shape determined from the 7r° and A 
sidebands. 335 ± 22 E+ events are obtained from the fit. 
We do a similar analysis to measure J/tp — > S^tt+A. The 
signal for E~ and the fitting result are shown in Fig. 5 
(b). The fit yields 254 ± 19 events. 

The efficiencies for J/4> — > E+7r~A and J/ip E^tt+A 
are determined to be 2.3% and 1.8% using 2 x 10^ 
MC simulated signal events, respectively. The branch- 
ing fractions are calculated to be B{J/^p E+tt^A) = 
(7.70 ±0.51) X 10-4 and B{J/^p E-tt+A) = (7.47 ± 
0.56) X 10^*, where the errors are statistical. The to- 
tal branching fraction of the two conjugate modes is 
B{J/tp -> E+TT- A -I- c.c.) = (15.17 ± 0.76) x lO""*. 



Background determination and upper limit on the 
number of signal events 
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FIG. 5: (a) M{7v°p) of J/i^ — > 7r*^p7r A candidate events from 
the signal region of Fig. 4 and (b) M{-k'^p) oi J/tp Tr^pyr^A 
candidate events. Dots with error bars are data, the solid 
histograms are the best fits described in the text, and the 
dashed histograms are backgrounds estimated from A and tt" 
sidebands. 



3. Measurement of J /ip E^vr A -)- c.c. 



The J/i! -> E^ 



^ K + c.c. events, where E"*" — s- Tr^p and 



A TT+p, have the same final states as the signal channel 
AAtt". Candidate events are required to satisfy <15, 
in addition to the common selection criteria in Section III, 
except for the A and A mass requirements. Figure 3 is a 
scatter plot of M(7r+p) versus M{tt~p) invariant mass for 
data and MC simulation. The two bands are the J/tp ^ 
E+TT" A -I- c.c. events. In order to select J/ip ^ E+tt-A, 
M{T:-p) > 1.15 GeV/c^ is required. Fi gure 4 shows the 
scatter plot of M{'j'-f) versus M(7r+p). The intersection 
region (central box) of the tt" and A bands corresponds 
to the J/tp ^ E+TT^A signal. The dots with error bars 
in Fig. 5 (a) show the distribution of M{tt'^p) invariant 
mass of the events in the central box (|M(77) — M(7r'^)| < 



Using the branching fractions for J/tp — > E+tt A + c.c. 
measured above and branching fractions available in the 
PDG [9], we obtain 29.2, 14.3, 14.2, 125.0, and 11.9 
background events from J/tp S°S", J/t/j E°E°, 
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FIG. 6: Invariant mass distribution of M(77) for J/tp 
AA77 candidates (dots with error bars) and normalized back- 
grounds (solid histogram). The dashed curves shows the tt" 
signal from MC simulated J/tp ^ AAtt''. The arrows denote 
the region of the TfO Si gnal defined in the text. We use different 
histogram styles to indicate leading backgrounds from J/tp ^ 
E°7r°A (-I-CX.) (circles), J/tp S+vf-A (+c.c.) (squares), 
J/tP E°E° (trmngles), J/tp E(1385)"S(1385)" (stars) 
and J/tp E°E° (rhombi), which contribute 46.0 ± 5.4, 
11.2 ± 1.3, 1.9 ± 0.4, 1.0 ± 0.3, and events in the defined 
tt" region. 
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J/jp S(138_5)0E(1385)°, J/^ T,°Tr°A (+c.c.), and 
J/V' S+7r~A (+C.C.) for the J/i/) AAtt" selection, re- 
spectively. We also studied backgrounds from other pos- 
sible channels listed in the PDG [9] that might contami- 
nate the TT*^ signal, but their contamination was found to 
be negligible. The histogram in Fig. 6 shows normalized 
backgrounds from all background channels. The normal- 
ized M(77) distribution of the background events is in 
reasonable agreement with the data. The dashed line in 
the figure shows the 7r° signal from MC simulated J/4> ^ 
AAtt^. To estimate the expected number of signal events, 
we define the tt" mass region as |M(77) — Af (7r°)| < 0.045 
GeV/c^, which is indicated in the figure and selects most 
of the tt" signal events. The numbers of tt*^ events in the 
mass region are found to be 54.0 ± 7.4 and 60.1 ± 9.5 for 
data and normalized backgrounds, respectively. By using 
the POLE method [13, 14], the upper limit on the num- 
ber of tt" events from J/-0 — > AAtt'^ is calculated to be 
11.2 at the 90% confidence level (C.L.). 



B. J/V? — »■ AAry 

Candidate events with two or three good photons arc 
selected, and the xlc required to be less than 15. Since 
the momenta of A and A are low in this channel, no 
requirement is made on the decay lengths of A and A; 
otherwise the efficiency would be extremely low. This 
is demonstrated in Fig. 6, where a decay length require- 
ment is made and no rj signal is seen. Figure 7 shows 
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FIG. 7: The AA invariant mass distribution for J/i/) ^ AA77 
candidates (dots with error bars), background from MC sim- 
ulated J/tp ~^ E^E" (hatched histogram), background from 
MC simulated J/Vj ^ E°E° (dashed histogram), and MC 
simulated J/ip AArj signal (solid histogram). The back- 
grounds are normalized according to the branching fractions in 
the PDG and the AAr; signal is normalized using the branch- 
ing fraction measured in this paper. The arrow indicates the 
M(AA) requirement, and events below the arrow are selected 
as J/tp —> AArj candidates. 



the AA invariant mass distribution after the above selec- 
tion. To remove the backgrounds from J/-0 — > S°E° and 
J/tP ^ S^SO, events with M(AA) > 2.55 GeV/c^ are 
rejected, since for the signal process, they are kinemati- 
cally prohibited . Dots with error bars in Fig. 8 show the 
invariant mass of A/ (77), and a clear rj signal is observed. 

To investigate possible backgrounds, we consider the 
following channels with A or S production: J/ip ^ 7AA, 
SOsO, S(1385)0S(1385)0, E^E°, 2(1530)05°, Y.°n°A + 
c.c, and S"'"7r~A -I- c.c. Using available branching frac- 
tions of these decay modes, we obtain 7.8, 27.6, 6.2, 
and 20.4 background events from J/i/j S°!E°, J/ip — > 
S(1385)0E(1385)0, J/V' ^ 2(1530)05°, and EOttOA-K c.c. 
in the mass region A^(77) > 0.4 GeV/c^, respectively. 
The background contribution from the AA sidebands 
(|M(7rp) - 1.141| < 0.01 GeV/c^) is evaluated to be 
3±2 events. Contamination from other possible channels 
listed in the PDG [9] that might contaminate the 77 signal 
is negligible. The shaded histogram in Fig. 8 shows the 
normalized backgrounds from the above channels. We fit 
the 77 invariant mass distribution with a MC simulated 
signal shape and a second order polynomial background. 
The fit yields 44± 10 events with a statistical significance 
of 4.8cr. 



C. ip(2S) AAtt" and AAr] 

The selection criteria for these two decays are similar to 
those for J/ip decays. A 4C kinematic fit to the hypoth- 
esis ipi'^S) ~* 777r"^7r^pp for candidate events with two 
good photons is performed, and the xlc required to be 
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FIG. 8: Fit to the 77 invariant mass distribution of J/tp 
AA77 candidate events selected in Fig. 7. Dots with error bars 
are data, the hatched histogram is the normalized background 
from all the channels considered, and the solid histogram is 
the fit to data using a histogram of the signal shape from MC 
simulation plus a second order polynomial for background. 
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less than 15. Backgrounds from ip{2S) tt+tt J/ip are 
rejected with the reqmrement |M^f "l' - M(J/'(/')| > 0.04 
GeV/c^, where M^^™!.' is the recoihng mass of tt+tt". 
Figure 9 (d) depicts the invariant mass distribution of 
the charged tracks. The peak around 3.1 GeV/c^ is 
from ip{2S) neutral + J/tp. In order to veto such 
background, |M(AA) - M{J/iP)\ > 0.05 GeV/c^ is re- 
quired. Furthermore, to suppress the background from 
'0(25') — > shown in Fig. 9, the invariant mass of AA 

is required to be less than 3.3 GeV/c^. 

Figure 10 shows the 77 invariant mass distribution af- 
ter the above selection, and we see no significant tt" or 77 
signals. In order to estimate the number of signal events, 
we define the signal regions as 0.09 < Af(77) < 0.18 and 
0.50 < M(77) < 0.60 (GeV/c^) for 7r° and rj, respectively. 
The number of signal events is found to be 4 in both re- 
gions. To estimate the backgrounds from the sidebands 
of n° and r?, (0.03 - 0.08) and (0.19 - 0.25) (GeV/c^) 
are taken as the sidebands of 7r°, and (0.43 - 0.49), (0.61 
- 0.67) (GeV/c^) are taken as the sidebands of 77. The 
numbers of background events from the sidebands of tt^ 
and r] are estimated to be 1 and 1.5. With the POLE 
method [13, 14], the upper limits on the numbers of sig- 
nal events at the 90% C.L. are calculated to be 7.0 and 
7.6, respectively. 



V. SYSTEMATIC ERRORS 



The systematic errors on the branching fractions are 
mainly from the efficiency differences between data and 
MC simulation in the MDC tracking, particle identifica- 
tion (PID), photon detection, kinematic fitting, the A ver- 
tex finding, and the decay length requirement and the un- 
certainties on the total number of J/iJj and ipi'^S) events. 
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FIG. 9: A/(AA) distribution for MC simulated events: (a) 
ipi2S) AAtv", (b) tp{2S) AA?7, and (c) ip{2S) 
(d) The M(AA) distribution for data. The arrow denotes the 
selection M(AA) < 3.3 GeV/c^ 
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FIG. 10: The 77 invariant mass distribution for candidate 
tp{2S) 77 AA events. Dots with error bars are data, and 
the histograms are MC simulated signal events. The arrows 
indicate the signal region of tt*^ and r; described in the text. 



The MDC tracking and particle identification (PID) 
systematic errors are estimated from the difference of 
the selection efficiencies of protons and antiprotons be- 
tween data and MC simulation [15]. The efficiencies 
are measured using samples of J/ip 7r+7r~pp and 
V'(2S') 7r+7r~pp, which are selected using PID for three 
tracks, allowing one proton or antiproton at a time to be 
missing in the fit [15]. The efficiency difference between 
data and MC simulation for one proton is from 2% to 5% 
depending on the proton momentum of the decay chan- 
nels. 



The photon detection efficiency is studied using J/^ 
P^tt" in Ref. [16]. The results indicate that the systematic 
error is about 2% for each photon. Therefore, 4% is taken 
as the systematic error on the photon efficiency for all the 
decays. 

The uncertainty due to the kinematic fit is studied us- 
ing many channels which can be selected purely without 
a kinematic fit [15-17]. It is found that the MC simulates 
the kinematic fit efficiency at the 5% level for almost all 
channels tested. Therefore, we take 5% as the systematic 
error due to the kinematic fit. 

The A reconstruction systematic errors are studied us- 
ing J/tp ^ AA [7, 8]. The A secondary vertex finding 
gives a systematic error of 0.7% for each A vertex, and 
the decay length requirement contributes 1.7%. The to- 
tal percentage error arising from A and A vertex require- 
ments is 3.7%. 

The systematic error of the background jihape can be 
determined by fitting the observed E+, and 77 sig- 
nal events with different background shapes. For J/tp — > 
I]+7r~A and J/tp li~7r+A, the background shape in 
fitting the E+ and S" is changed to a second order poly- 
nomial. The differences in the numbers of fitted S"*" and 
S" events are found to be 2.2% and 1.5%, respectively. 
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TABLE I: Summary of systematic errors (%). 



Source 


J/-!/) ^ AAtt" 


J/V" ^ AAri 


i){2S) AAtt" 


tp{2S) AAtj 


J/iP 


>-A J/ip E-TT+A 


Tracking and PID 


7.0 


14.0 


6.0 


12.0 


7.0 


6.0 


Photon efficiency 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Kinematic fit 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


A vertex requirement 


3.7 












Background shape 




3.0 






2.2 


1.5 


Number of good photons 


3.0 


3.0 


3.0 


3.0 






Total number of events 


4.7 


4.7 


4.0 


4.0 


4.7 


4.7 


Total 


11.6 


16.6 


10.1 


14.5 


10.8 


10.1 



TABLE IL Measured branching fractions or upper limits at 90% confidence level (C.L.) for all the studied channels. Here, 
B{A ir~p) = 63.9%, B(E+ ^ 7T°p) = 51.6% and Bij) 77) = 39.4% are taken from the PDG. 



Channels 


Number of events 


MC efiiciency(%) 


Branching fraction ( X 10 ^) 


J/iP AAtt" 


< 11.2 


0.75 


< 0.64 


J/%j) AAri 


44± 10 


1.8 


2.62 ± 0.60 ± 0.44 


ip{2S) AAtt" 


< 7.0 


2.5 


< 0.49 


ip{2S) AAri 


< 7.6 


2.9 


< 1.2 


J/iP E+TT-A 


335 ± 22 


2.3 


7.70 ±0.51 ±0.83 


J/V) E-TT+A 


254 ± 19 


1.8 


7.47 ±0.56 ±0.76 


J/V) E+7r"A + c.c. 






15.17 ±0.76 ±1.59 



For J/ip AAry, the background shape is ehanged from 
a second order polynomial to a first order one, and the 
difference in the number of fitted signal events is about 
3%. 

The uncertainty caused by the requirement of two good 
photons is estimated by considering the percentage of 
events without fake photons in the sample oi J/ip ^ AA. 
It is found that the difference in the percentages of events 
without fake photons between data and MC simulation 
is 3%, which is taken as the systematic error for the re- 
quirement of two good photons. 

Finally, the results reported here are based on a total 
of 58 M J/tp events and 14 M tjj{2S) events. The un- 
certainties on the number of J/tp and ■ip{2S) events are 
4.7% and 4.0%, respectively. Table I lists the systematic 
errors from all sources. Adding all errors in quadrature, 
the total percentage errors range from 10% to 17% for all 
the studied decay channels. 



VI. RESULTS AND DISCUSSION 

Table II lists the results for J/tp and tp{2S) decay into 
AAtt" and AA77, as well as J/ip E+7r-A + c.c We also 
list the total branching fraction for the conjugate modes, 
where the common systematic errors have been taken out. 



Except for J/ip AAtt" and J/ip S+tt-A + c.c, the 
results are first measurements. Interestingly, the result 
of J/tp — > AAtt^ presented here is much smaller than 
those of DM2 and BESI [3, 4]. In previous experiments, 
the large contaminations from J/tp — > S'^tt'^A + c.c. and 
J/tp S+TT-A ± c.c. were not considered, resulting in a 
large value of branching fraction for J/tp AAtt^. The 
small branching fraction of J/tp — > AAtt" and relatively 
large branching fraction of J/tp AAry measured here 
indicate that the isospin violating decay in J/ tp decays is 
suppressed while isospin conserving decays are favored, 
which is consistent with expectation. 
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